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ABSTRACT 

We address the dynamical evolution of an isolated self-gravitating system with two 
stellar mass groups. We vary the individual ratio of the heavy to light bodies, jj, from 
1.25 to 50 and alter also the fraction of the total heavy mass Mh from 5% to 40% of 
the whole cluster mass. Clean-cut properties of the cluster dynamics are examined, like 
core collapse, the evolution of the central potential, as well as escapers. We present in 
this work collisional iV-body simulations, using the high-order integrator NBODY6++ 
with up to A/"* = 2 • 10 4 particles improving the statistical significancy of the lower- 
A/** simulations by ensemble averages. Equipartition slows down the gravothermal 
contraction of the core slightly. Beyond a critical value of \x m 2, no equipartition 
can be achieved between the different masses; the heavy component decouples and 
collapses. For the first time the critical boundary between Spitzcr-stable and -unstable 
systems is demonstrated in direct TV-body models. We also present measurements of 
the Coulomb logarithm and discuss the relative importance of evaporation and ejection 
of escapers. 
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1 INTRODUCTION 

The internal evolution of a star cluster in dynamical equi- 
librium is governed by its tendency to achieve a thermal 
velocity distribution (i.e., Maxwellian with energy equipar- 
tition) through small changes of velocity during 2-body en- 
counters between stars, a phenomenon dubbed relaxation. 
The relaxation time is the average time after which a star's 
moving direction has been defle cted by 90 degrees relative to 
its original orbit JSpitzerlll987lh Relaxation produces major 
changes in the structure of the cluster while keeping it in 
dynamical equilibrium. 

The collapse of the central core makes up an important 
phase and probably the most fascinating aspect of the dy- 
namical cluster evolution. There are three mechanisms act- 
ing in different ways in order to achieve the collapse: equipar- 
tition, evaporation and gravothermal instability. In a real 
cluster of stars all of these processes are present, but in the 
idealized models, the ones we shall describe, it is possible to 
isolate the specific processes and gain some understanding 
of t he particular e ffects. 

ISpitzerl dl969T) set about an analysis on segregation of 
masses in globular clusters systems that would lead later 
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to a broad ensemble of different analyses and techniques. 
For some clusters it seemed impossible to find a configu- 
ration in which they have dynamical and thermal equilib- 
rium altogether. The heavy component sink into the centre 
because they cede kinetic energy to the light one on the 
road to equipartition. In most of the cases, equipartition 
happens to be impossible, because the subsystem of mas- 
sive objects becomes self-gravitating before. Thermal en- 
ergy flows from the inner part to the outer regions. Whereas 
the outer regions of the cluster do not alter significantly 
their temperature, the inner regions, the core, loses heat 
and, so, contracts and becomes hotter. A self gravitating 
system has a negative thermic capacity. This phenomenon 
has been obser ved in a big number of works using dif- 
ferent methods jHcrion 1973 : lHenonlll975l: ISpitzer fc Shut 
19751: ICohnlll98ol: iMarchant fc Shapirolll980l: IStodolki 
ITakahashil f 



19821 : iTakahashil 119931: iGiersz fc Heari? 
1995HSpurzem fc Aarsethll996l:|MakinJl996l:lQuinlanll99el : 
Drukier et al.ll999l:lJoshi et alj|200(l etch The late phase of 
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core collapse is the same as for a single mass model, because 
the heavy components do not interact with other stars any- 
more. 

There is an a mple evidence for mass - segreg ation 
in observed clusters. | Mc Cmujhr ean fc Staufferl i 1994) and 
iHillenbrand fc HartmanrT " i^)9«) provided a new deep in- 
frared observations of the Trapezium cluster in Orion that 
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clearly show the mass segregation in the system, with the 
highest mass stars segregated into the centre of the clus- 
ter. This is a clear-cut evidence for mass-segregation of stars 
more massive than 5 Mq toward the cluster centre and some 
evidence for general mass segregatio n persisting down to 1- 
2 in the Orion Nebula cluster. iRaboud fc Mermilliodl 
( 1998) study the radial structure of Praesepe and of the very 
young open cluster NGC 6231. There they find evidence for 
mass segregation among the cluster members and between 
binaries and single stars. They put it down to the greater 
average mass of the multiple systems. 

At this point, the question looms up whether for very 
young clusters mass segregation is due to relaxation, like in 
our models, or rather reflects the fact that massive stars are 
formed preferentially towards the centre of the cluster, as 
some models predict. 

To answer such questions there is a clear necessity for 
models that give us an accurate description of the evolution 
of multi-mass models based on direct-summation numerical 
schemes. ^^^^^^^^ 

The simplest case of a bimodal mass spectrum is a start- 
ing point to take care of. This is a relatively good app roxi- 
matio n if stellar black holes are the heavy component llLed 
Il995l) . Such two-mass simulations are exclusively studied in 
this work. 

ISpitzerl (1 19691) gave an analytical criterion to determine 
whether a two-component system can, in principle, achieve 
energy equipartition or not. According to his analysis, en- 
ergy equipartition between the light and heavy component 
can exist if 



S := 



3/2 



< 0.16 



(1) 



Where Mi and Mh are the total masses in light and heavy 
stars and mi and mh their individual masses, respectively. 
Spitzer's work was based on many strong simplifying as- 
sumptions. 

A number of authors has addressed the problem of 
thermal and dynamical equilibrium in star clusters from 
a nu merical point of view, with direc t A-body simula- 
tions JPorteeies Zwart fc, McM illan 200(J, Monte Carlo sim- 
ulationT lSpitzer^^Iar^^9^i r a nd with direct integration 
of the Fokker-P lanck equation Jlnaeaki fc Wivantd Il98l 
iKim et al]|l998f) . As regards the Monte Car lo scheme, re- 
cent and very detailed numerical calculations (matters et alJ 
I2000T) have suggested a different criterion, 



( Mh 
\Mi 



mh 
mi 



< 0.32 



(2) 



The limitations inherent in this approach motivate us 
to embark on more accurate models of this scenario with 
the help of A-body methods, where Newtonian gravity is 
essentially treated without approximation s. 

In his pioneering work, Ivon Hoernerl l|l960F ) performed 
calculations with A/"* = 16 particles on the best comput- 
ers available at that time. Rapid improvements in computer 
technology (both hard- and software) facilitated larger as 
well as more accurate calculations. The amount of 10 4 par- 
ticles was reached by dSpurzem fc Aarsetbligg'ri) . and paral- 
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Table 1. Overview of the simulations. A model is described by 
q and fi and assigned to a capital Roman letter. RND = random 
setup; INS = all heavy masses placed inside; OUT = all heavy 
masses outside. As regards II, was inside and ouside. The 
series VI with A/"* = 5 ■ 10 3 can be envisaged as a model for the 
Orion Nebula 



lei machines and special purpose computers do even manage 
fifty times more nowadays. 

In this article we describe the simulation models based 
on these methods as well as their initial conditions. We 
present the results of a wide parameter space, which has 
been explored by direct A-body modelling for the first time. 
We also extract the important parameters describing the 
core collapse and equipartition of energies. 



2 ORGANIZATION OF THE SIMULATIONS: 
NOMENCLATURE 

For all simulations in this work we employed a Plummer 
sphere model in global virial equilibrium. The particles are 
treated as point masses without softening of the gravita- 
tional force, but with regularisation of close encounters in- 
stead. We only consider two different mass species as the 
most simple approximation of a realistic mass spectrum. 
Since this analysis aims to isolate the essential physical pro- 
cess of mass segregation, we ignore stellar evolution, clus- 
ter rotation and a tidal field as well as primordial binaries; 
binary formation occurs only during the late stage of the 
evolution and do not effect our objectives. 

In our notation, a model will be determined by its frac- 
tion of the heavy mass component, q := Mh/M c i, and the 
mass ratio of the individual particles, n := m^/mi. The 
model is assigned to a capital Roman letter. Each model 
consists of a number of runs that differ only in their random 
number seed which produces different initial setup of posi- 
tions and velocities of the particles for the same distribution 
function. The runs are physically equivalent. Models making 
up a logical unit for comparison are gathered to a series, see 
Table □ 

After the work bv llnagaki fc Wivantd dl984f) . we con- 
centrate on q = 0.1 (Series I), since this is the value which 
had the fastest evolution, as they proved, and study the 
evolution for a wide range of /i's (from 1 to 50). This choice 
is guided by observations in youngest star clusters like the 
Trapezium in the Orion N ebula Cloud, where a mass rang e 
of w 0.1-50 M is found Jffillenbrand fc Hartmanr3ll998li . 
This cluster also exhibits clear mass segregation that cannot 
be explained by the simple theory of "general" mass segre- 
gation driven by two-body relaxation. Given the extreme 
youth of the stars with highest mass, a primordial segre- 
gation has been suggested, in wh ich they were formed a t 
locations close to the dense centre feonnell fc Davies!ll998f) . 
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We examine the time scales for a random setup of parti- 
cles (RND) by assigning a mass mi and mh to each body. For 
(j is a relative quantity, we fix m\ to unity, and vary mh in 
steps that are given in the top row of Table H The two mass 
populations have A/i and Ah members (we will elaborate on 
this ahead), whose spatial coordinates whose positions and 
velociti es are picked up ran domly according to a Plummer 
sphere jAarseth et al.lll974l) . 

In Series II, we place the heavy particles completely 
either in the centre (INS) or in the outskirts (OUT) and 
compare the evolutionary patterns with the random setup of 
Series I. In Series III, IV and V we alter the fraction of heavy 
masses, q, to test for this parameter. The full parameter 
space is graphically illustrated in Figure Additionaly, we 
performed simulations of a special configuration, which is 
related to a mass ratio in the Orion Nebula Cloud (Series 
VI). 

In order to reduce the statistical noise, we per- 
formed a large number of runs and averaged the data 
set into an ensemble-model. The statistical quality of 
such an ensemble-averaging is comparable with one sin- 
gle calculation containi ng the full particle set of all runs 
jGiersz fc Heggiejl994aTl . This means that, since the fiducial 
model is a cluster with A/"* = 2.5 • 10 3 particles, we carried 
out 20 runs and the whole set containing Af* = 5 • 10 4 has as 
little noise as a high-A/"* model. In Series I, we also performed 
50 runs for A/"* = 10 3 , and 10 runs for A"* = 5 • 10 3 . Addi- 
tional models containing 10 4 particles (4 runs) and 2 ■ 10 4 
particles (1 run) were also performed. 



2.1 Parameters of two— mass models 

In order to describe the physics of our models, we have to 
define a set of three parameters: 



RND: Models A-G with N = a500 



A"* = Ai + Ah, 



/i = mh/mi, 



q = M b /Md. (3) 



Since the number of the force calculations per crossing time 
scales with A/"* , computational efforts restrict the choice of 
A/"*. We can express A/i as follows, 

Al h m h A"h 
q M\+ Alh ~~ miA/i + m h A"h 

mi/iAh 
m\N\ + mi/x/Vh 
MA"* - A/i) 
Ai + M (A; - Ai) ' 

The absolute number of light particles is then 



(4) 



Ai = 



(1 - q)nM<. 



(5) 



q - m + m 

Some authors (e.g. Inag aki fc Wivantolll984ll define a 
slightly different parameter q := Alh/AIi. With this defini- 
ton, equation turns out to be A/i = uAf* Kg + p,) ■ In the 
follow ing, we employ the notation of ISpurzem fc TakahashJ 
Jl995T) , as defined in (J^J . The absolute numbers of the heavy 
particles, Ah = A/"* — A/i, are shown in Table[3]for the models 
of Series I (q = 0.1). /i is given in brackets on top of each 
column. 

A different way of fixing the mass ratio /i is by means 
of the average mass of the stars: 




each model: aven-uie of 20 Runs 



Figure 2. Comparison of the shrinking core radius in the models 
A— G of Series I given in Nbody units (see Appendix [EJ 



where (m) — A4 c i/A/"* and rm is the i-th component in a 
multi-mass cluster. With g; = A4i/A4 c i, we have for the 
general case of k different mass components 



q»A/"* 

fj-i 



i)qiAf* 
rrn 



Mi 
mi 



(7) 



particles in the i-th mass bin. 

The advantage of this expression lies in the simpler han- 
dling, if more than two masses are present. In the equation 
0, (k — 1)! parameters of rrn/mj would be necessary for k 
mass components, while employing (m) reduces the amount 
of the a priori definitions of m^'s to k — 1. 



3 RESULTS FOR AI h /Af c 



0.1 



fit = rrii/{m), 



(6) 



In this section we deal with the situation in which the frac- 
tion of heavy stars makes up 10% of the whole cluster mass. 
We investigate the physical processes occurring from a ran- 
dom initial distribution and compare them with previous 
literature on this subject. The essential facts about the evo- 
lution are qualitatively visible in Figure [5] The core radius 
shrinks with time, and the cluster collapses under its self- 
gravity. The time scale for the collapse is the shorter the 
larger the mass ratio /i between heavy and light stars is. 
The formation of binaries in the core stops the collapse and 
allows its re-expansion. We focus our attention on the vari- 
ations of the evolutionary processes for different values of 



3.1 Data of ensemble averages 

We vary fi from 1.25 to 50.0 and obtain the core collapse 
time, t C c, for each run by considering two values: first, the 
moment of the minimum core radius, t (r c | m i n ) and, second, 
the deepest central potential, t ($| m m). The output data was 
written for each A-body time unit. Because of large fluctu- 
ations between two subsequent data points, we applied a 
"sliding average" over r c and $, based on the following al- 
gorithm 
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Figure 1. Parameter space of the mod- 
els examined in this work. The position 
of the Roman letters indicate a model 
determined by its q and fi. The dotted 
line is the boundary of the equip arti- 
tion stable region after SDitzer Jl969l) . the 
dashc d-dotted line is th e stabil ity criterion 
after iLightman fc Fall ll97Sf> . the dot- 
dot— dot— dashedline isthe prediction by 
llnagaki Wivant 3 <1984h and the dashed 
lin e is the empirically p roposed condition 
bv lWatters et alJ fcOOOft 
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Table 2. Absolute numbers of heavy stars for models of Series I. Next to each model we give the value for /i in brackets 



Evolution for various N 




Figure 3. Evolution of the central potential of selected runs of 
Model A for A/"* = 10 3 , 2.5 ■ 10 3 , 5 ■ 10 3 , 10 4 and 2 ■ 10 4 particles 



^ W — l 

TZi = - V A 

111 — ^ 



— w/2 i 



(8) 



where TZi is the resulting value of the original variable At 
averaged among w neighbouring data points. The choice of 
w — 5 turned out to be best fitting one. 

The average of the times found from the minima of r c 



and <3> is defined to be the core collapse time of the run: 

t Tc . + tif> 



(9) 



Figure 0] illustrates the distribution of the data points 
of r c and $ m i n at their corresponding core collapse times for 
four selected models. The logarithmical time-axis was cho- 
sen to show the relative scatter between the point clouds. 
When looking at a particular run, the core collapse times are 
likely to show a significant discrepancy between the determi- 
nation from the potential or from the core, i.e. t$ does not 
necessarily correspond to t Ic , but the averages of both, (t Io ) 
and {£*), yields a good concordance for the whole model. 
Thus, (t cc ) is a good value to characterize the core collapse 
time of the model, 



(tec) = 



J2[(tr c +U)/2] _ {t rc ) + {U) 



(10) 



uns ^ 

In this last equation A/runs is the number of runs and we 
obtain an equality between the left-hand and right-hand side 
by linearity. 

As we can see in Fig. (|1J, the relative variance of the core 
collapse times increases for smaller values of (t cc ) . The mean 
values for the eight models of our Series I are summarized in 
Appendix A. The errors given there are the standard devia- 
tion from the runs' mean t cc , divided by the square root of 
the number of the runs "vWruns- The error of (t cc } is roughly 
2-5% for most of the models, consistent with the relative 
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Figure 4. Models A, C, D, and F of Series I with JV* = 2500. 
Each point represents the minimum value of the central potential 
<f> c (upper panel) and core radius r c (lower panel) at the time 
when these parameters attain their minimum. The average values 
for all models are listed in the Appendix 



errors determined by Spurzem & Aarseth (1996) and the 
half-mass evaporation times by Baumgardt (2001). 



4 EVOLUTION OF THE CORE RADIUS 

The evolution of the core radii for the models A-G is plotted 
m Figure Each curve is an ensemble average of 20 runs. 
For models with /j, approaching unity, the core radius shrinks 
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up 
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0.7855 
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0.5035 


0.4325 



Table 3. Measured exponent for core radius expansion in the 
post-collapse phase. See text for further explanation. 



as in an equal-mass system (at late collapse times, to the 
right). A linear time scale (not shown here) suggests that the 
collapse phase sets in when the core radius has contracted 
to about 25% of its initial radiu s . This is in agreement with 
the results of lGiersz fc Heggid (Il994al) . For high /i's, this 
happens from the start of the simulations. The rapid con- 
traction of r c is due to the very massive stars falling quickly 
to the centre. The contraction stops at higher r c -values than 
for the low-/i models, and a quick expansion of the core fol- 
lows. 

The behaviour at the moment of core bounce is illus- 
trated in Figure^] where the minima of the potential and 
the core radius are plotted versus [i. For a fixed A/* and 
values of /i between 3-10, the core collapse is carried out 
by approximately the same number of particles, but they do 
not draw together as close as for the equal-mass case or very 
small /i's. Therefore, the density and the central potential 
are less deep than for /i < 2, and the minimum of the core 
radius is not so profound. At high /i's, the effect is reversed: 
Even some few heavy components are massive enough to 
deepen the potential. The collapse itself is less distinct, as 
seen from the shallower r c in the right panel. In this range 
of values for /t, the depth of the potential results from the 
combination of these two effects: Heavy masses build up a 
strong gravitational field, but their kinetic motion does not 
allow a long-lasting vicinity. Equipartition becomes then im- 
possible. The differences of the depths for various N* are 
discussed in section [H] 

For a system in which the central object has a small 
mass and the energy production is confined to a small central 
volume we have that the core-radius r c should fol low an 
expansion proportional to the power-law of t and t cc < HenorJ 



expansion proportional to trie power-law ol t and t cc I 
19651: IShapirdtl977l iMcMillan et"ailll98ll: iGoodmanl 



mi), 



OC (t — tcc) 



2/3 



(11) 



We measured the slopes of the expanding r c 's in Figure 
|5|by fitting two straight lines embracing the fluctuating data 
of each model to construct an upper and a lower margin. 
These lines give two independent measurements for the slope 
and are presented in Table [3] and they show a decreasing 
trend for higher /i. 

The lower lines seem to be steeper than the upper ones 
for the most of the models. On the other hand, the simu- 
lations reach different stages of the post-collapse evolution, 
and the expanding branches exhibit different lengths. Espe- 
cially, the low-/x models (A and B) are not far advanced for 
precise measurements, while the high-p models (F and G) 
appear distorted about the time of collapse such that the 
ons et of the self-similar ex pansion is difficult to find (see 
also lGiersz fc Heggielll994lJl . From the theoretical point of 
view, there is no argument for a different behaviour of the 
core expansion when unequal masses are present. 
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Figure 5. Minima of potential (upper panel) and core radius 
(lower panel) at the moment of core bounce. The collapse is car- 
ried out by a less number of core particles when fi rises. The filled 
circles correspond to the single-mass case. See text for explana- 
tions in more detail 



5 CORE COLLAPSE TIMES 

As we explained in the introduction, two-body relaxation 
causes star clusters to redistribute the thermal energy 
among stars. Since this kind of heat transfer acts on the 
relaxation time scale, a core collapse is similarly ensued in 
gravitationally unstable systems. 

The core collapse time is best studied numerically. For 
equal-mass models it ranges about t cc m 330t rc (0) during 
the sel f-similar phase , or about 12-19 half-mass relaxation 
times (lQuinlanlll996l) . In this article, Quinlan gives a time 
scale of 15.7 t r h for an isolated cluster, if an isotropic veloc- 



ity distribution is assumed. iTakahashj |l995l) also modelled 
Plummer spheres but for an anisotropic case, and deter- 
mined the collapse time to about 17.6 t T h- Other authors 
find similar factors, and we shall adopt 



17.5 tr 



(12) 



Though the values are used in most studies of the core 
collapse, they are a poor guide for clusters with a mass spec- 
trum, e.g. globular clusters have central relaxation times 
that are typically ten, sometimes a h undred times s horter 
than their half-mass relaxation times fouinlanlll996t ?). 

So far, the core collapse time is only found empirically 
from a large number of numerical simulations, for there ex- 
ists no analytical theory which would predict it a priori from 
cluster properties, e.g. the star number, IMF, or concentra- 
tion parameters. 

It has been suggested that the nucleus consisting of two 
mass-components collapses in a ti me shorter than the eq ual- 
mass cluster by a factor of 1/fi (iFregeau et alJ 120021) . In 
broad terms this is so because the equipartition time is short- 
ened according to the time-scale for equipartition iSpitzerl 
Il969l) . The upper panel of Figure |S| shows detailed calcula- 
tions in a wide /x-rangc of this. The mean times, {t cc ), have 
been plotted versus p for the complete sample of our models 
in Series I. 

In the upper panel, the core collapse times converge 
smoothly to the collapse time for an equal-mass cluster, 
tcc.i, as n approaches unity. Between /i ss 2 and 10, a some- 
what linear decline is visible, but far beyond the stability 
boundary, run towards a constant value, since binaries are 
likely to have more importance the larger /x is. In the same 
Figure, in the lower panel, we display the factor of "col- 
lapse acceleration" for various /Li's. This gives us a mea- 
sure to which percentage a cluster with two masses evolves 
faster than the single-mass case. The dashed line is the fi -1 - 
decline, i.e. a cluster with two mass-species having a ratio 
H would collapse /x-times earlier than its single-mass equiv- 
alent according to 



1 

^cc. (A, OC tcc,l- 



(13) 



At fi = 2, the decoupling of the equipartition-based 
instability and the gravothermal instability seems to take 
place. For jj, — > 1, the evolution occurrs more slowly 
than the l//i-decrease, because the tendency to equipar- 
tition drives the initial evolution and slows down a "purely" 
gravothermal collapse. This seems to suggest that as the 
mass difference between mh and m\ becomes less impor- 
tant, the system collapses like a single-mass cluster. Be- 
yond the critical value [i = 2, an early decoupling of the 
two mass populations occurs. The heavy components try 
to reduce their large velocity dispersion, but they rapidly 
accumulate in the centre and interact preferably with them- 
selves. As a consequence, equipartition is harder and harder 
to achieve and the evolution proceeds only due to the re- 
distribut ion of heat within the two ; almost separated com- 
ponents toettwieser fc In agaki 198J3). The light stars evapo- 
rate out of the core and take away the thermal energy to the 
outskirts, while the heavy components increase their binding 
energy. The latter ones collapse like a single-mass subclus- 
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Figure 6. Upper panel: Core collapse times of the models in 
Series I. The times shorten non— linearly with increasing ji. The 
dotted l ine at fi = 1.275 is the equipartition criterion given by 
ISpitzerl Il969t Eq. 3-54) for a value of Xmax = 0.16 and q = 0.1. 
The core collapse times of the corresponding equal-mass clusters 
are indicated by filled dots. The error bars are smaller than the 
symbols. In the lower panel all collapse times are normalised to 
tcc,l (the core collapse time for the single— mass case) The simula- 
tions for N-t, = 2.5 ■ 10 3 are connected with a solid line for clarity. 
The circles at left are the collapse times for an equal— mass cluster 
according to eq. 1121 . and the dotted line is the stability boundary 
by Spitzer with Xmax = 0.16 <Spitzerlll987l Eq. 3-55). 



ter. The release of energy transported away and the heat 
transfer works effectively and leads to an accelerated col- 
lapse. 

If /x 3> 1, the situation turns into a case of dynamical 
friction: A significant fraction of particles are drowned into 
a homogeneous sea of light stars and, like in an ordinary 
frictional drag, their motion suffers a deceleration. It is in- 
structive to see that only for large A/"*, the absolute number 
of heavy stars seems sufficient to maintain the linear slope of 



accelerated cluster evolution. The slope follows the dashed 
line a bit longer before bending towards that constant value. 



6 MASS SEGREGATION 

The process of mass segregation for the six models A-F is 
illustrated in the Figure where the mean mass of the stars 
that are inside a specified Lagrangian shell, is shown (i.e., 
a Lagrangian shell is the volume between two Lagrangian 
radii, which contai n a fixed mass fraction o f the bound stars 
in the system; see lOiersz fc HeggieHl994al) . With the light 
and heavy masses randomly distributed, each shell exhibits 
the same average mass at the beginning. In the course of the 
evolution, the inner shells assemble the heavy bodies, and 
raise the mean mass. The half-mass radius and the outer 
shells lose their heavy stars rather quickly and remain be- 
low the value for the average mass, because the light ones 
outnumber the heavy components significantly. The cluster 
is stratified by mass. 

The segregation of masses pr opagates simultaneously 
with the contraction of the core. iGiersz fc Heggiel dl996l) 
have already noticed the self-similarity of this process: In 
the models A and B we can clearly see how the inner layers 
decouple stepwise one after the other until the final stage 
is reached. After the collapse, the profile does not change 
much; the effect has only influence on the early evolution- 
ary phase, not in the post-collapse. The moderate decline in 
the models E and F can be explained by the heavy particles 
escaping: Because of their relatively small absolute number 
(see Table |5J, they get ejected from the core and leave the 
cluster so that the surplus of small bodies depletes the mean 
value in the shells. A fraction of high-mass stars in the cen- 
tre interact strongly in a few-body process and kick out each 
other, and the core is gradually "evacuated" from the heavy 
stars. In our data we find an enhanced fraction of high-mass 
cscapers occuring immediately after the collapse confirming 
this scenario; many of them exhibit enormous kinetic en- 
ergies and some even escape as bound binaries (see Sec. [TJ . 
However, some heavy latecomers enter the inner regions, but 
a balance is found between the incoming and outgoing mass 
flux. In general, the dynamical processes in the core do not 
influence the properties of the cluster as a whole. 

The mean mass in the 1%-Lagrangian shell never at- 
tains its full "capacity", i.e. the value of that could prin- 
cipally be gained if this shell was completely populated by 
heavy stars; which is only possible if the heavy stars are 
51%. It means that there is always a number of low-mass 
stars entering and leaving the innermost region that keep 
(m)i% at a constant fraction of the maximally attainable 
level. This level is at about 0.8 if fi is small, and drops 
to 0.5 (Bh for the highest /i's. 

Another illustration of the segregation process is given 
in Figure |H] It shows the shrinking distance moduli of 18 
heavy-mass stars of one typical run of a model with a high 
mass-ratio (Model G, = 25, Run no. 4, A/"* = 5 ■ 10 3 ); this 
model consists of 22 heavy mass stars initially, but 4 stars 
escaped at some time before the moment for which we show 
the plot and they were excluded from it. 17 dashed lines were 
overlapped demonstrating how the orbits of most particles 
draw rapidly closer. The core collapse occurs at tec — 33.5. 
A solid line lifts off one particle of example that was knocked 
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Figure 7. Average mass in Lagrangian shells for models A— F show the stratification of masses in the cluster. The segregation of heavy 
masses proceeds in agreement with the global evolution of the cluster. The mean mass is indicated 
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G05_RND/Run04: Segregation of heavy particles 
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Nbody time 



5 ■ 10 . Each of the particles has a mass 
_ 25 times larger than a light particle (or- 
100 bits not shown). The core collapse of this 
run occurred at t cc =33.5 



out of the chaotic region due to a close encounter, but its 
kinetic energy was not sufficient enough to leave the system 
and it segregated inwards again. One other particle, whose 
pericenter was at 0.9 iV-body-radii, remained in the halo 
and did not take part in the collapse, but it sank to inter- 
mediate distances at about t cc ~ 70. 



7 ESCAPERS 

The estimation of escape rates is often based on idealized 
models, whose simplifying treatments sometimes lead to dif- 
ferent results. Th e comp lexity of this topic is reviewed by 
iMevlan fc Heggiel (^997). For the results that we show in 
this work, we should mention that: 

• Theories based on diffusive or small-angle relaxation 
phenomena yield a different escape rate than theories in- 
volving individual two-body encounters; the former is often 
denoted as "evaporation", the latter is related to "ejections", 
which dominate for isolated clusters, like here 

• The rate of escape is not a constant, while the evolution 
of the system proceeds, even in the pre-collapse phase 

• An increasing concentration in the core as well as the 
growth of anisotropy tends to enhance the escape rate 

• Furthermore, the escape rate is strongly mass depen- 
dent; different mass spectra and segregation alter it 

• A tidal field lowers the energy threshold for escape 

• A sufficient abundance of binaries (both, primordial and 
formed) has a substantial effect on high-velocity escapers 
which take away energy from the system 

In view of these complications, we have to take care in the 
interpretation of the data. We show in this section TV-body 
simulations investigating the rates for various particle num- 
bers and the variation on /i, in particular. 

An escape is defined by a particle having both positive 
energy and its distance from the density centre exceeding a 
limiting radius. We have chosen the distance to be 20 times 
the half-mass radius, r^. 



The particles were removed from the calculation when 
both conditions were fulfilled; we shall call them "removed 
escapers". We focus here on escapers occuring before the 
time of core collapse, t cc . One run contains a number of 
particles which have gained positive energy but not reached 
the distance for a removal yet. Such kind of particles domi- 
nate when t cc is very short. In particular, the high-^i mod- 
els collapse within a few tens of A-body-time units, and a 
large number of particles, that are going to escape, would 
be missed. This s ituation resembles the "energy cut-off" 
jBaumgardtll200ll) . and we shall call these particles "po- 
tential escapers". Whether some of them will be scattered 
back to become bound members again or really escape, is 
a complex process that is out of the scope of our subject. 
iBaumgardtl i200lT) estimated that a fraction of 2% of the 
potential escapers might return to the system. As a first ap- 
proach to the general properties of escapers, we will assign 
to the number A4sc all removed escapers plus potential es- 
capers at the epoch of t cc . Because of the scatter of the core 
collapse times (section l3.11 . each run was checked separately 
for its escapers that occured before that run's individual core 
collapse time — not the ensemble's (t cc ). (7V CS c) is the aver- 
age among the individual Apse's. 

Figure |5] gives a typical example for the energy distri- 
bution of all removed escapers in one run of Model C with 
A/"* = 5 ■ 10 3 ; the potential escapers are disregarded here. 
The energy is plotted against the time of removal; it is mea- 
sured in units of kT — moQ av , where m is the average mass 
and ao,av the average ID velocity dispersion in the core, 
measured for the initial model. Light bodies are represented 
by filled dots, heavy ones of = 2.0 by diamonds. The 
collapse occured at t cc = 643.5 for this run. In the begin- 
ning, the light particles diffuse slowly from the system with 
small energies (evaporational effect). After the core collapse, 
the mass dependence is more complicated, for a second class 
of escapers joins: high-energy particles, whose energies are 
higher by two orders of magnitude. A fair fraction of the es- 
capers are heavy stars. As in the statistics for equal-masses 
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bv lGiersz fc Heggid dl994aTl . it is natural to associate them 
with ejected stars that go back to three-body interactions 
in the very centre. Mass segregation has widely finished at 
the time of core collapse, and interactions in the core start 
depleting the high-mass population. Other runs of the same 
model exhibit a similar picture of the physical scenario. 

The analysis of the models in regard to fi reveals in- 
teresting views on the escape mechanism. A summary is 
shown in Figure [TU1 and the legend for the symbols is given 
in panel (e). Panel (a) gives the fraction of escaped stars, 
(A/'esc}/A/'*. In a single-mass cluster (/i = 1), about 2.5% of 
the stars leave the system before it collapses. When intro- 
ducing a second mass, this fraction drops to 0.2-0.5% until 
H w 3. The reason is that the rate of escape (panel e) is 
nearly constant for small /i, but the shorter collapse times 
cause a smaller progress of the escape mechanism, and thus 
a smaller A/" CS c- When massive bodies mi, 3> rn\ are present, 
a larger fraction of stars receives positive energy and turns 
into potential escapers. 

Panel (b) gives the ratio of the removed escapers to 
the potential escapers as described above, (-/V rem ) / (A/pot) • 
For n close to the equal-mass case, the core collapse time is 
large. Thus, the accumulated number of removed particles is 
larger by « 10 times than the number of potential escapers 
at the moment of t cc . For high n's, removed particles are 
scarce, while potential escapers did not have time to cross 
the cluster and turn into removed ones. Therefore, the po- 
tential escapers make the overwhelming majority. The slope 
of the decrease is — 1 , and it is similar to the shortening of 
t cc when /j, increases. The vertical dependence on A/+ mir- 
rors the increasing amount of removed escapers due to the 
longer t cc -times for larger A/"*. 

Panels (c) and (d) deal with the individual masses of the 
escapers. When one heavy star escapes, it raises the average 
escaper mass. In the models fi > 10, the absolute number 
of high-mass escapers is zero for most of the runs. This is 
the converse situation of panel (a): For massive stars, the 
escape is difficult. Just 1 or 2 heavy components appear and 
take away a significant fraction of mass from the cluster. 
They are rather evaporated objects than ejected in a close 
encounter, for their energies are relatively small. Usually, the 
heavy escapers occur in the late post-collapse phase and do 
not receive our attention here. From panel (d) and fi ^ 5, we 
find that fa 3% of all escapers are heavy components which 
leave the system before t cc . 

Panel (e) presents the rates of escape within a crossing 
time for different /i. The computation was adopted from 
Wielen (1975) as 



dt 



= U 



(Ncsc) 



(14) 



where the brackets denote averages among the runs 
within a model and t CT is the average number of escapers 
per crossing time. For equal-mass clusters our rate is ~ 
0.2 — 0.3 in accordance w ith other A-body results listed by 
iGiersz fc Heggid lll994aD . The slight decline until /i ~ 3 and 
the stronger increase afterwards is in exce l lent a greement 
with the theoretical expectation bv iHenonl Jl969T) (Fig. 1). 
Though his models differ from ours in q, the branch of the 
curve related to our models exhibits the same shape. We 
interpret the curve such that two different mechanisms pro- 



duce escapers in the pre-collapse phase under consideration: 
Relaxation dominates if the individual masses do not differ 
much, thus evaporation causes a steady mass loss from the 
system. When fi is increased slightly, the effect of evapora- 
tion cannot advance so far. The escape rate is reduced then. 
For higher /i's another mechanism takes over: Massive stars 
exhibit a strong gravitational focusing. The gravity of one 
single heavy particle attracts more light stars and its energy 
is distributed among them. So, a multitude of light stars eas- 
ily gains positive energy and heads for escape. The frequent 
two-body encounters of one heavy particle leads rather to 
ejections involving high energies than an accumulation of 
small escape energies. 

Panel (f) of Figure [TU1 shows the mean energy carried 
away by the stars. The near constancy suggests an indepen- 
dence on fi, but a subdivision into removed and potential es- 
capers (not shown here) reveals differences between the two 
groups. For removed escapers, the mean energy increases as 
much as a factor of 10 over the whole p-range. This confirms 
the ejection scenario explained for panel (e). On the other 
hand, the potential escapers show a constant but somewhat 
lower energies on average. Since they make up the larger 
fraction, the mean is depleted in the high-^i regime. 

Finally, we can conclude that models resembling the 
equal-mass model tend to lose their mass by a slow evapo- 
ration process, while energetic ejections outweight in high-^i 
models. At n ~ 3, both processes are exchanging their dom- 
inant role. After mass segregation has come to a stop near 
core bounce, three-body encounters in the core start deplet- 
ing the population of m,h-stars. 



8 LARGER PARTICLE NUMBERS 

Some dependencies on different particle numbers have al- 
ready been tacled in the foregoing sections. Here we present 
a direct comparison of the cluster evolution for the A-models 
first. This model is close to the uniform mass case and has 
the longest evolution time. Figure |3] shows the raw data 
of the minimum potential for five runs, i.e. we applied no 
smoothing or averaging. 

The most obvious feature, is the increasing core collapse 
time in pretty accordance to the increasing relaxation time 
proportional to A/*/ln(7A/+) (see section 8.1). Note that the 
fluctuations of the data in the pre-collapse phase are smaller 
for higher A/"*, because the global potential is smoother. In 
the post -collapse, the fluctuations are nearly the same, for 
the number of core particles, N c , is of the same order for 
each of the five runs. 

The assumed constancy of A/" c (where the subscript "c" 
stands for core) leads to a second topic concerning the dif- 
ferent amplitudes of the potential minima. This has already 
been mentioned in connection with Figure |3 though the 
upper panel of Figure 1111 points to the variations on Af+ in 
a more clear way. For models resembling the equal-mass 
case, the maximum depth is a function of A/"*. The reason 
is that the collapse is only halted when the rate of energy 
production in the core becomes the same as the energy rate 
going out via th e heat flux of the gravothermal instability 
llGoodmanlll987ft . The outflowing energy is produced by the 
formation of binaries in three-body encounters, and this be- 
come important when the density is sufficiently high. As the 
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Figure 9. Model C (fj. = 2.0), Run no. 04, 
with A/"* = 5 ■ 10 3 . Energy of an escaping 
particle is given in units of kT and plotted 
2.0 2.5 against its escape time scaled to the core 

collapse. Diamonds denote heavy masses 



core radius of a large-TV* cluster contains a larger number 
of stars initially, it has to get rid of almost all of them. The 
final Nc at core bounce is nearly a constant and the evolu- 
tion of the core radius advances deeper in order to provide 
the den sity necessary for binary formation. iGiersz fc Heggiel 
described this scenario in terms of the fraction of 
core radius to half-mass radius, r c /rh'- the larger A/*, the 
smaller that fraction. This is visible for models with /x close 
to unity in Figure ^3 Though, the lower panel suggests a 
slight increase of core particles for large N models, the ab- 
solute values of A/" c stay in the order between 10 and 40. 

The situation looks somewhat different for models with 
high fi, in particular for fj, = 25 and 50. The central po- 
tential shows its deepest value when some few of the heavy 
masses have gathered in the core, while the light stars have 
not changed their density distribution. The presence of the 
massive stars causes the deeper potential, and it appears not 
as profound as for nearly equal masses (see Figure|5]for com- 
parison). Therefore, the amplitude of the collapse, ("i&mm}, 
depends on A/" c only weakly. 



8.1 The Coulomb logarithm 

Now, we will focus on the coefficient 7 in the Coulomb log- 
arithm, which is of relevant importance for the evolution 
of the star cluster, for it is intimately connected with the 
time-scale associated with th e evolution of the cluster, the 
relaxation time (Eq. 2-62 of Spitzer 1987). The value of 
this quantity, In (7A/"»), has been estimated to be 7 = 0.11 
iGiersz fc Heggiall994af) for a single-mass cluster, but the 
variations on \i are not known precisely. One possible way 
to determine this quantity is by comparing the evolution of 
the same model but for different values of A/"*. 

The upper panel of Figure 1121 shows the ensemble- 
averaged Lagrangian radii of our Model A for two differ- 
ent A/"*. For the time t, the value of a Lagrangian radius in 
the 2500-body model was determined, and then the corre- 



sponding time t' at which the same value was reached in the 
5000-body model. 

The ratio of these two times is a scale factor, Sf, that 
should be also equal to the ratio of the relaxation times for 
the different A/"* before core rebounce: 



(15) 



i_ = s A^ln( 7 A^) 
t : ~ A^ln( 7 A^)' 

By repeting this procedure for each time step and for each 
Lagrangian radius we get an Sf that is plotted in the lower 
panel of Figure H2l Before the core collapse, the time ratios 
show a remarkable constancy (disregarding the initial set- 
tling period), and the inner Lagrangian radii have a good 
agreement. The scale factor in the constant pre-collapse 
range is Sf = 1.675 for this model. With the definition 
v := A^/A/"*, a re-arrangement of eq. (1151 yields 



7 = 



AT* 



V /{S t -v) 



(16) 



We computed the 7's for each model using the data 
sets of A/"* = 2.5 • 10 3 and A/* = 5 • 10 3 . The results are 
given in Tableland plotted in Figure The first two rows 
of the Table define the /^-model, the third gives the scale 
factor as determined from the figures analog to !12l the fourth 
is the error measured from the widths of the fluctuating 
lines in the stable regime (horizontal part), the fifth is the 
outer exponent in Eq. 1161 with v = 2 and the sixth are the 
resulting 7's. 

The comparison of the data sets "2.5- 10 3 — > 5T0 3 " pro- 
vides a more accurate scale factor than the other data sets 
with lower A/"*, because they rest upon a higher statistical 
significancy: Our set with A/"* = 10 4 can be averaged over 4 
runs only, and the 10 3 -set is biased to low-A/"* physics. 

Very small variations in Sf cause large fluctuations in 
the two exponents of eq. 1 Kit and alter 7 significantly; the 
formula is very sensitive to both, Sf and v. The error was 
computed as the difference of an "upper" and a "lower" 7 
that results from the thickness of the Sf-line. 

As for different v, we present another example illus- 
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Figure 10. Summary of escaper data for the different models, (a) Fractions of removed escapers before the core collapse plus potential 
escapers at t cc . (b) Ratio of removed escapers to potential escapers. (c) Mean mass of escapers. (d) Relative number of m h -escapers 
among all escaping particles, (e) Escape rates, (f) Mean energy of escapers in kT. The error bars are usually omitted, though they are 
given for the JV* = 2.5 ■ 10 3 set (triangles) except in the panels (c) and (d); the error bars are almost invisible in (c), and larger than the 
panel size in (d), because the number of high-mass escapers varies a lot among one model. 
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model A B C D E F G 

H 1.25 1.5 2.0 3.0 5.0 10.0 25.0 



St 


1.675 


1.800 


1.785 


1.750 


1.775 


1.845 


1.680 


ASf 


0.045 


0.010 


0.035 


0.030 


0.175 


0.045 


0.190 


u/(S { - u) 


-6.154 


-10.00 


-9.302 


-8.000 


-8.889 


-12.903 


-6.250 


7 


0.0142 


0.2048 


0.1263 


0.0512 


0.0948 


1.5321 


0.0152 



Table 4. Scale factors and j'e for the models with different fi. The data were determined by comparing the data sets A/"* = 2500 and 
Ml = 5000, giving a particle ratio v = 2 



trating the big variations when determining the Coulomb^y. 
Table |^| shows the results for the equal-mass models using 
various particle numbers A/"* and A/"* . The Sf shows a reason- 
ab le behaviour for d i fferent v's and is conform to the result 
bv lGiersz fc Heggid <ll994al) shown in the last column. The 
resulting 7, however, differs at least by a factor of 5. 

The value 7 = 0.11 mentioned above is actually found 
due to the comparison of A r -b ody and Fokker-Planck simu- 
lations of equal-mass models llGiersz fc Heggielll994al . Fig- 
ure 6). From their A-body-A-body simulations with parti- 
cle numbers of A/"* = 500 and Ml = 2000 (their Figure 5), it 
is easy to extract an Sf = 2.95 and get a 7 — 0.098. For the 
same particle ratio, v = 4.0, we obtain a 7 that is half of 
theirs (last two columns in Table But our absolute par- 
ticle numbers are 5 times larger. We can conclude from our 
analysis that the choice of M* and Ml essentially contributes 
to the Coulomb-7. 

Finally, we are left with a fair range of possible values for 
7. An estimate indicates that 7 ranges somewhere between 
0.01 and 0.20. A similar range was given by Giersz & Heggie 
(1996) for the case of a power-law mass spectrum: 0.016 < 
7 < 0.26. 



9 VARIATIONS OF THE MASS FRACTION 

So far we have presented results on mass segregation for 
various values of /1. This section deals with three addi- 
tional series, in which we alter the fraction of the heavy 
masses, q. Such ki nd of study has been discussed by 
Watters et al 1 fcooot with a Monte Carlo approach or by 
InageJti^^W^antol (1198^1 . who simulated two-component 
clusters by means of Fokker-Planck modelling. They fixed 
H to 2.0 and investigated the core collapse times as well 
as the achievement of equipartition. They showed that the 
evolution of the central potential in units of the half-mass 
relaxation time, t T h, for 4 different q's is fastest when q ~ 0.1. 

We present in this section the question whether equipar- 
tition can be achieved between the two components while 
segregation is on work and a comparison of the TV-body data 
with the above mentioned literature. 

We investigate the equipartition parameter 



Series q /^crit 



mi a f 



(17) 



which gives the ratio of the kinetic energies between the 
heavy and light stars in the core. At the start its value is 
about jj, and heads for unity. When both mass species find 
a state of energy equipartition, £ = 1 is reached, and we call 



III 0.05 2.49 



0.10 2.03 



IV 0.20 1.87 

V 0.40 1.75 



2.032 


(M+ 


= 1000) 


2.048 


(A4 


= 2500) 


2.197 


(M* 


= 5000) 


2.111 


(AT* 


= 10k) 


1.762 


(AT, 


= 20k) 



Table 6. Values of fi at which £ m i n > 1.05. Equipartition cannot 
be attained anymore. See text for further details 



the system equipartition stable ( after ISpitzerilT969Tl . other- 
wise a £ m in indicates the closest approach to it. 

When examing a particular run, the data of £ is very 
noisy, especially, for small-A/"* simulations. The main cause 
for this is the small number of particles inside the core ra- 
dius. For this reason, we decided to evaluate for the particles 
inside twice the core radius, and then apply the smoothing 
procedure of section l3.ll We define £ m i n = 1, if eq. 1171 drops 
below unity at any time during the evolution; otherwise we 
set it to the deepest peak. Figure fTTI shows an example for 
the parameter £ from two models with M* = 104. The B- 
model does reach equipartition at t = 1157, so £ m i n is set 
to 1; the F-model approaches to it down to £ m i n = 2.95 at 
t = 106. 

The average from all the runs' particular £ m i n 's is de- 
noted by (£min) and taken as the most reliable data for the 
model. The (£ m in)'s are plotted versus n in Figure IT51 The 
symbols belong all to the Series I with q = 0.1, but a different 
particle number. The solid line connects the M± = 2.5 • 10 3 
values (triangles). The additional Series III (q = 0.05), IV 
(q = 0.2), and V (q = 0.4) are shown as solid lines. Error 
bars are given only for the four sets with M* = 2.5- 10 3 ; they 
are smaller than the symbols in most cases. 

As expected, the graphic shows that equipartition takes 
place for small p's, but when \i becomes significantly greater 
than 2, (£ m m) recedes from unity. At /J, = 2, about half 
of the individual runs with q = 0.1 did succeed to reach 
Cmin = 1, at least for a moment. Those runs, which did not 
find a state of full equipartition, tried to reduce the kinetic 
difference halfways to the core collapse, but then departed 
shortly after the closest approach. 

That figure provides a good basis for the judgement on 
equipartition stability. It is obvious that (f m i n ) varies for dif- 
ferent fractions of heavy masses, q: The less amount of heavy 
masses in a cluster, the closer equiparti tion is reached. This 
is consistent with the results by both llnagaki fc Wivantol 
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Eq. Ms 



1000 -> 2500 1000 -» 5000 2500 -s- 5000 2500 10,000 Giersz fc Heggie (1994a) 



<5f 

- i/) 

7 



2.5 
1.96 
-4.630 
0.028 



5.0 

3.21 
-2.793 
0.018 



2.0 
1.60 
-5.00 
0.006 



4.0 
3.09 
-4.396 
0.044 



4.0 
2.95 

-3.810 
0.098 



Table 5. Comparison of our equal— mass models with four different particle numbers and the iV-body-analysis by lOiersz fc Heggiel 
Jl994al) in the last column 



Closest approach to equipartition 



10 - 



l - 




10 



For Series I: 

x N = 20,000 

<> N=10,000 

□ N= 5000 

A N= 2500 

* N= 1000 



100 



Figure 15. The averaged minimum of 
the equipartition parameter £ for stars in- 
side twice the core radius. When both 
mass components are in equipartition, 
Smin equals 1. The solid lines connect the 
data points for q = 0.4 (Series V), q = 0.2 
(Series IV), q = 0.1 (Series I with trian- 
gles), and q = 0.05 (Series III). These four 
simulation series were made with A/"* = 
2.5 • 10 3 particles, and error bars are given 
for them. The dotted line is an arbitrary 
threshold for equipartition stability at 1.05 
as explained in the text 



l)l984h and lWatters et alJ l)2000ri . The latter explore an even 
wider range to very low g's down to 0.0015 (their set "B"). 

In order to check the stability border, we look now for 
a "critical /x" at which equipartition stability is not given. 

Instead of assuming that equipartition is happens for 
(Cmin) = 1, we will give a small tolerance for this and define 
the poi nt of "turning away" from equipartition at (£mm) = 
1.05, as lWatters et all <200d) did. It is denoted by the dotted 
horizontal line in Figure Hoi By linear interpolation between 
the lower-next and upper-next data point, one obtains that 
this threshold is exceeded at the points fj, CI i t given in Table 

m 

The value for Series I was averaged from all simu- 
lation sets for different A/"*. All (£ m in)'s are plotted Fig- 
ure 1161 and fitted for a direct compar ison of the Monte 
Carlo results from IWatters et all d200(j) . The positions of 
our /i cr it are marked b y a small filled dot. Th e apparent 
differ ence for the line bv | Light,man fc FalJ l|l978T) is just ba- 
cause IWatters et al defined their fraction of heavy 

masses as q = Mh/Mi, while we use q = A4h/A4 c i (see sec- 
tion . The graphics and the results are consistent with 
each other and ease the comparison. 

Although we have four points for /i C rit in order to check 
the equ ipartition bounda r y, the y follow the theoretical func- 
tion bv lLightman fc FalJ l|l978F) in fair agreement. 

The formula suggested bv lWatters et all l)2000f) cannot 
be ruled out, for it is based on models in a low-q regime, 
which is difficult to access with our A^body simulations. 



Ho wever, our exp erimental data to show that the criterion 
by ISpitzerJ il969l) for the stability boundary, appears too 
strong, especially at mass ratios /i close to 1. 

/icrit ~ 2 is the point which was already recognized 
as the transition of an equipartition-dominated and a 
gravothermal-dominated collapse (section Below /i C rit, 
the core collapse proceeds slowlier than the theory pre- 
dicts, because equipartition governs the initial phase. When 
/i > Merit , the gravothermal instability always wins the com- 
petition between the two effects: Before the thermal equilib- 
rium can fully be achieved, the massive stars have already 
segregated to the centre and collapse independently from the 
light stars. With the total fraction q being high, the massive 
component appears almost self-gravitating and is decoupled 
from the beginning. 

iBettwieser fc Inagakil lll985l) showed that there is an 
"optimal" q that favours the collapse of a two-component 
system. If q is small, the collapse of the heavier component 
proceeds in the external field of the lighter stars, and the 
collapse time decreases rapidly. On the other hand, if q is 
very large, the amount of thermal energy contained in the 
core of the massive stars is so large that the core of the light 
ones cannot receive all the redundant energy. Hence the ten- 
dency to collapse is postponed, because the heat gained has 
to be dumped away into the halo and this p rocess is sl ower 
than the tendency to equipartition; see also llLeelll995l) . Fi- 
nally, the extreme cases q — > and q — > 1 approach the 
single-mass value for a negligible contribution of the heavy 
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Figure 11. Upper panel: The mean potential minima at the mo- 
ment of core bounce as as function of the total particle number 
A/"*. Models close to the equal-mass case show a deeper peak when 
A/* is rising than the high-/x models. Explanation is given in the 
text. Lower panel: The mean core number increases only slightly 
with A/"*. Disregarding the two highest /n, the number of minimum 
core particles is between 10 and 40 



and light component, respectively. This is the reason why 
the core collapse time attains a minimum at a certain frac- 
tion of the heavier stars (Bettwieser & Inagaki 1985). The 
JV-body models confirm this (Figure [T71 , 



10 CONCLUSIONS 

We have studied mass segregation and equipartition in ide- 
alized star clusters containing two mass components. The 
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Figure 12. Upper panel: Five selected Lagrangian radii of Model 
A for A4 = 2.5 • 10 3 and JV* = 5 ■ 10 3 . Lower panel: Scale factor Sf 
for this model computed from the comparison of the Lagrangian 
radii 



results provide potential insight into the evolution of young 
star clusters, dynamics of star forming regions, and the de- 
gree of initial segregation which is needed to match the ob- 
servations. Also, this study is of theoretical interest to check 
the classical ideas on equipartition and the thermodynamic 
behaviour of self-gravitating systems. 

The parameter space was widely analysed for both pa- 
rameters using Plummer models: the ratio of individual star 
masses, y, — m^/mi, that was varied from 1.25 to 50; and 
the fraction of the total heavy masses, q = Mh/M c i, that 
was chosen to be 0.05, 0.1, 0.2, or 0.4. We also directed 
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Figure 13. Upper panel: Scale factors measured from figures 
analog to !12l for various fi. Lower panel: Coulomb-7 as determined 
from JV-body simulations with different particle numbers from 
cq. 1161 . The error bars are larger than the panel size 
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Figure 14. Equipartition parameter for two runs from model B 
(upper panel) and F (lower panel), respectively, containing 10,000 
particles. The data was smoothed out with a smooth width w = 5 
as described in Chapter 13.11 



our attention to the valu e q = 0.1, which was proved by 
llnagaki fc Wivantol lll984l) to exhibit the fastest evolution 
time. For this fraction, we modelled different particle num- 
bers, with A/"* ranging from 10 3 to 2 • 10 4 . For all other cases 
we used a particle number of TV* = 2.5 • 10 3 . The statisti- 
cal significancy of the data was considerably improved by 
ensemble-averages: A large number of physically equivalent 
runs differing just by a random seed for an initial setup was 
performed, and the data base was gathered for an overall 
average. 

For the first time we have empirically checked the bor- 



ders between an eq uipartition stable and unstable configura- 
tion JSpitzerll969l) with the accurate direct iV-body method. 
His criterion involves the parameters q and fi and is given in 
Eq. Q. Moreover, we determined the ranges of the longest 
and shortest core collapse time by comparing different initial 
setups of the particles. 

From the results of our simulations we draw the follow- 
ing conclusions on the global cluster evolution: 

(i) The evolution of a self-gravitating system with two 
mass-components depends strongly on the ratio of the in- 
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Equip art it ion — stability 



i ^ □ □ 




Figure 16. Closest approach to equipartition for the models in 
our parameter range of /i (mass ratio of heavy stars to light stars) 
and q (mass in heavy stars relative to total cluster mass, see 
also Figure^. The figure has been adjusted (symbols, axes, and 
lines are accommodated as explained in the text) s o that one can 
easily compare it with the results with those of Wattcrs ct al. 
fcOOOl) (Fig. 6). The solid line is their stability criterion, Eq. 
the straight dashe d line is f rom [Smtzeij |l 969l) . and the curved 
dashed line is from lLightman jTFalir^97«) 



^=2.0 



Figure 17. Core collapse time t cc as a function of q. The in- 
dividual mass ratio, /j, wa s fixed to 2.0. The gas models by 
iBettwieser &; Inagakil ll985l) show a broad minimum for the core 
collapse times in terms of the relaxation time. The value of q for 
fastest evolution is rj 0.1—0.2. The results of the TV-body mod- 
els for the four series are discussed in the text. The time on the 
vertical axis is given in TV-body units 



dividual masses, /i = rah /mi. For q = 0.1, when ji is larger 
than ~ 2, the heavy masses fall to the centre and reduce the 
relaxation time proportional to 1/fi. Subsidiary, the core col- 
lapse time is shortened in the same way (Fig.0. Smaller /x's 
go into equipartition which slows down the evolution. Large 
/j,'s give rise to a small number of heavy particles, and the 
situation turns into a process of dynamical friction. 

(ii) There is a competition between equipartition of en- 
ergies and gravothermal instability. If /j, is close to unity, 
equipartition hinders the gravothermal collapse in the ini- 



tial phase, but will never prevent it. As the heavy stars 
congregate in the central regions, they decouple from the 
light component and perform their own collapse. When [i 
exceeds a critical value of about 2, equipartition can never 
be achieved (Fig. I15H . 

(iii) The boundary between t he stable and u n stable 
regimes is close to the theory by iLightman fc Fall dl978l) 
when q is > 0.2. For q lower than « 0.1, other works on 
equipartition JWatters et alJ l200Ct) seem to show a better 
agreement, but it is difficult to assess their range with N- 
body simulations because the absolute number of very heavy 
particles is too low in our standard A/* (Fig.ll M. When q be - 
comes larger than 0.1, the criterion given bv lspftzeil (ll96STt 
appears too strong, and his theory fails: According to Eq. Q 
no equipartition should be possible at all, while our mod- 
els as well as the Monte Carlo simulations bv lWatters et all 
teOOOl) show the contrary. The cause is probably that Spitzer 
assumed a global equilibrium when deriving the analytical 
formulas. His basic ideas on the processes, however, are still 
in fair accordance with our results. 

(iv) Segregation proceeds on the relaxation time scale si- 
multaneously with the evolution of the cluster. The maxi- 
mum level of segregation is attained when the core has col- 
lapsed and begins to expand. Then, the mass shells in the 
cluster have also adjusted to a stable balance of an in-going 
and outgoing mass flux (Fig|7Jl. After the core collapse, the 
degree of stratification remains constant. 

(v) The escape rate shows two branches along the /x- 
axis indicating two different mechanisms being on work 
(Fig. \Wk). For mass ratios resembling the equal-mass case, 
the escape seems to be governed by evaporational effects 
in the pre-collapse phase. Since the core collapse time de- 
creases for a moderately rising fi, evaporation does not ad- 
vance far, such that the mass loss is also reduced. For high 
^t's (> 3), we believe that escapes are rather a matter of 
ejections; the very massive particles distribute their kinetic 
energy to a large number of stars, and the escape rate in- 
creases roughly with /i. The details of these processes, how- 
ever, need more analysis tracing the path of the escapers or 
a countercheck with other simulation methods. As for the 
equal-mass case, the escape rate is consistent with other 
A-body simulatio ns jWielenlll975l : I Giersz fc Heggielll994al : 
iBaumgardt et al.ll2002l) . 

(vi) Light masses play an important role in the heat trans- 
fer from the core to the halo. Even a small fraction of them 
moderates the heat flux between the central source and the 
outer sink. The "optimal" fraction q of heavy to light stars is 
about 0.15-0.20 in accord with previous results from gaseous 
models (Fig.tTTl. 

Our results gained with the A-body method confirm 
previous simulations with other techniques (gas models, 
Monte Carlo), but also reveal hidden effects like the small 
"deceleration" of the gravothermal collapse due to equipar- 
tition. Further astrophysical assumptions are likely to alter 
the results, e.g. the introduction of primordial binaries in 
various fractions, a tidal field that accelerates the collapse 
as well as the dissolution of clusters, or rotation. Multiple 
mass components or a continuous spectrum make the anal- 
ysis more complex but are important for the understanding 
of the observations. Our variation of fi can principally be 
generalized by merging into a multi-mass model, but care 
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has to be taken when defining the parameters. The new pa- 
rameters appearing for the models are to be checked, e.g. 
the slope of the initial mass function (if a power-law is as- 
sumed), and how the ratio of the highest mass to the mean 
mass, m max /(m), influences the segregation. 
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APPENDIX A: APPENDIX: DATA TABLES 

The following tables contain the mean values of the simula- 
tions of star clusters and are organised as follows. The title 
specifies the Series, q, and particle number. 
Column 1: model name and /j 
Column 2: (t cc ), mean core collapse time 
Column 3: ("^mm), depth of central potential 
Column 4: (r c ), mean core radius in the moment of core 
bounce 

Column 5: (Ac), number of particles in the core at the time 
of smallest core radius 

Column 6: (£ m m), mean of the closest approach to equipar- 
tition 

Column 7: (A/" OS c), number of escaped stars until core 
collapse time 

Column 8: {E esc ), mean energy of escapers 



APPENDIX B: 
TIMESCALES 



iV-BODY-UNITS AND 



Dimensionless units, so-called "JV-body units", were used 
throughout the calculations. They are obtained when set- 
ting the gravitational constant G and the initial total clus- 
ter mass A id equal to 1, and the initial total energ y 
E to -1/4 iHeggie fc Mathieul Il986l: lAarseth et alJ 119741) . 
Since the total energy E of the system is E = K + W 
with K = ^M c \{v 2 ) being the total kinetic energy and 
W = -(3tv/32)GM^/R the potential energy of the Plum- 
mer sphere, we find from the virial theorem that 



2 64 R 



1 |http : //www . hlrs ■ de/\ 



(Bl) 



R is a quantity which determines the length scale of a Plum- 
mer sphere. Using the specific definitions for G, M c \, and E 
above, this scaling radius becomes R = 37r/16 in dimension- 
less units. The half mass radius can easily be evaluated 
by the formula fe.g. ISpitzeHll987t) : 



M(r) = M c i 



3 /i? 3 



(l + r 2 /7? 2 ) 3/2 



(B2) 



when setting M(r h ) = \M C \. It yields r h = (2 2/3 - 
l)~ 1/2 i? = 1.30 7?. The half- mass radius is located at the 
scale length of R — 0.766, or about 3/4 "A-body-radii" . 
The initial half-mass crossing time of a particle is 



GM 



5/2 



(2/5)3/2 ■ 

Since the A-body-unit of time, £nb , is 1 when 



£nb = 



GM 



;/2 



(-4£)3/ 2 ■■ 



(B3) 



(B4) 



immediately follows that t cr /fNB = 2 a/2- 

In the situations considered here, the evolution of the 
cluster is driven by 2-body relaxation. Therefore, a natu- 
ral time sca le is the half-m ass relaxation time. We use the 
definition of ISpitzerl lll987t) . 



0.138^ / R\ 



In A 



1/2 



GM C 



1/2 



(B5) 



For a Plummer model, the half- mass radius is R1/2 
1.305 R. Md is the total stellar mass. 
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Equal-mass: AC = 1000 

model (t cc ) {*mi n > (r c ) (AC) (£ min ) (ACsc) (E CBC ) 

EQMASS01 345.8 ±5.2 -4.72 ±0.10 0.0208 ±0.0020 11.2 ±0.4 26.6 ±1.1 0.88 ±0.18 



Equal-mass: AC = 2500 



model 


(tcc) 


(* min > 


(rc) 


<AQ 


(£min) (ACsc) {^csc > 


EQMASS02 


716.0 ±12.0 


-6.02 ±0.16 


0.0086 ±0.0006 


14.8 ±0.8 


64.2 ±3.1 0.89 ±0.24 



Equal-mass: AC = 5000 

model (t cc ) (*min> {To) (AC) (£min} (ACsc) (Ecsc) 

EQMASS05 1210.2 ±10.6 -7.58 ±0.27 0.0045 ±0.0003 17.7 ±0.8 117.1 ±5.0 0.71 ±0.26 



Equal-mass: AC = 10,000 



model 


(*cc) 


(*min> 


(ro) 


(AC) 


fenin) (ACsc) (E esc ) 


EQMASS10 


2312.9 ±34.7 


-8.75 ±0.47 


0.0025 ±0.0004 


21.3 ±2.2 


256.3 ±7.1 2.55 ±1.43 



Series I: q 


= 0.1, AC = 


1000 












model (fi) 


<*cc) 


(*min) 


(rc) 


(AC) 




(ACsc) 


(-Eesc) 


A (1.25) 


309.1 ±5.8 


-4.60 ±0.11 


0.0206 ±0.0009 


11.1 ±0.4 


1.002 ±0.002 


22.9 ±1.2 


0.70 ±0.17 


B (1.5) 


256.0 ±5.1 


-4.10 ±0.07 


0.0275 ±0.0011 


11.7 ±0.3 


1.011 ±0.005 


18.4 ±1.1 


0.70 ±0.11 


C (2.0) 


162.6 ±3.5 


-3.52 ±0.04 


0.0400 ±0.0016 


12.0 ±0.5 


1.043 ±0.010 


9.5 ±0.5 


0.78 ±0.24 


D (3.0) 


86.1 ±2.1 


-3.27 ±0.06 


0.0596 ±0.0026 


12.0 ±0.6 


1.264 ±0.033 


4.8 ±0.4 


1.02 ±0.26 


E (5.0) 


45.8 ±1.4 


-3.24 ±0.06 


0.0797 ±0.0040 


12.5 ±0.8 


1.892 ±0.071 


0.0 ±0.0 


0.00 ±0.00 


F (10.0) 


25.4 ±0.9 


-3.65 ±0.08 


0.0954 ±0.0039 


12.1 ±0.7 


3.426 ±0.187 


4.1 ±0.4 


1.28 ±0.33 


G (25.0) 


18.0 ±0.7 


-4.72 ±0.12 


0.0785 ±0.0030 


7.3 ±0.4 


4.992 ±0.237 


12.9 ±0.2 


1.09 ±0.12 



Series I: q 


= 0.1, AC = 


2500 












model (fi) 


<*cc) 


(*min) 


(rc) 


(AC) 


<?min) 


(ACsc) 


(-Eesc) 


A (1.25) 


677.8 ±8.1 


-5.51 ±0.12 


0.0100 ±0.0005 


15.9 ±0.7 


1.000 ±0.000 


61.3 ±2.0 


0.77 ±0.27 


B (1.5) 


558.4 ±8.8 


-5.02 ±0.11 


0.0137 ±0.0006 


15.9 ±0.8 


1.001 ±0.000 


45.2 ±2.2 


0.55 ±0.14 


C (2.0) 


366.6 ±8.4 


-4.06 ±0.11 


0.0207 ±0.0011 


14.4 ±1.0 


1.043 ±0.028 


24.9 ±1.6 


1.18 ±0.28 


D (3.0) 


180.1 ±2.9 


-3.50 ±0.10 


0.0336 ±0.0019 


16.2 ±2.3 


1.188 ±0.017 


9.7 ±0.6 


0.70 ±0.12 


E (5.0) 


86.1 ±3.3 


-3.17 ±0.10 


0.0567 ±0.0043 


17.0 ±1.4 


1.807 ±0.056 


7.3 ±0.8 


1.43 ±0.27 


F (10.0) 


40.2 ±2.0 


-3.45 ±0.15 


0.0811 ±0.0054 


19.9 ±1.8 


3.218 ±0.086 


8.0 ±0.7 


0.89 ±0.12 


G (25.0) 


21.8 ±1.1 


-4.11 ±0.18 


0.0905 ±0.0048 


19.7 ±1.7 


7.150 ±0.282 


12.4 ±0.9 


0.86 ±0.15 


H (50.0) 


17.9 ±1.1 


-6.16 ±0.36 


0.0792 ±0.0039 


11.7 ±1.0 


11.245 ±0.556 


25.2 ±2.6 


0.98 ±0.15 



Series I: q = 0.1, AC = 5000 



model (fi) 


(tcc) 


(*min) 


(rc) 


(AC) 


«min> 


(ACsc) 


(Eesc) 


A (1.25) 


1183.7 ±17.7 


-6.92 ±0.18 


0.0055 ±0.0004 


19.8 ±0.7 


1.000 ±0.000 


103.9 ±5.0 


0.44 ±0.13 


B (1.5) 


993.8 ±10.3 


-5.85 ±0.12 


0.0071 ±0.0003 


16.9 ±1.0 


1.000 ±0.000 


76.9 ±3.2 


0.57 ±0.14 


C (2.0) 


642.6 ± 4.4 


-4.62 ±0.20 


0.0117 ±0.0011 


15.8 ±1.4 


1.016 ±0.008 


38.0 ±1.9 


0.47 ±0.09 


D (3.0) 


331.5 ±11.6 


-3.72 ±0.08 


0.0204 ±0.0013 


16.8 ±1.4 


1.189 ±0.037 


19.3 ±2.6 


1.08 ±0.45 


E (5.0) 


157.9 ± 6.3 


-3.36 ±0.12 


0.0302 ±0.0025 


14.9 ±1.1 


1.715 ±0.039 


14.2 ±1.9 


1.38 ±0.27 


F (10.0) 


64.3 ± 3.5 


-3.11 ±0.09 


0.0742 ±0.0055 


29.4 ±2.8 


3.738 ±0.256 


10.4 ±1.0 


2.17 ±1.01 


G (25.0) 


33.3 ± 0.9 


-4.43 ±0.34 


0.0842 ±0.0025 


29.1 ±1.4 


7.394 ±0.214 


17.7 ±1.2 


0.97 ±0.09 


H (50.0) 


22.4 ± 0.7 


-5.64 ±0.43 


0.0789 ±0.0089 


26.2 ±3.5 


12.962 ±0.358 


25.8 ±3.1 


1.15 ±0.18 
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Series I: q = 0.1, A/"* = 10,000 



model (fi) 


(tec) 




<* mi „> 


(rc) 




<€min) 


(A4 ac > 


(Eesc ) 


A (1.25) 


2169.9 ±24.7 


-8.35 ±0.52 


0.0027 ±0.0003 


20.0 ±2.7 


1.000 ±0.000 


227.3 ±3.9 


0.40 ±0.11 


B (1.5) 


1886.8 ± 


8.6 


-7.05 ±0.39 


0.0042 ±0.0001 


20.8 ±1.0 


1.008 ±0.006 


172.3 ±9.3 


0.35 ±0.06 


C (2.0) 


1218.0 ±20.1 


-5.85 ±0.37 


0.0053 ±0.0007 


15.8 ±1.8 


1.035 ±0.014 


91.3 ±4.3 


0.78 ±0.17 


D (3.0) 


595.1 ± 


8.5 


-3.97 ±0.17 


0.0107 ±0.0010 


15.0 ±1.9 


1.170 ±0.024 


35.3 ±3.4 


1.03 ±0.27 


E (5.0) 


277.6 ± 


9.7 


-4.18 ±0.32 


0.0235 ±0.0004 


19.8 ±1.4 


1.700 ±0.029 


19.3 ±1.4 


0.63 ±0.08 


F (10.0) 


123.5 ± 


2.1 


-3.80 ±0.31 


0.0351 ±0.0055 


22.3 ±3.9 


3.065 ±0.108 


25.5 ±2.0 


1.19 ±0.04 


G (25.0) 


52.4 ± 


2.3 


-4.92 ±0.42 


0.0741 ±0.0108 


36.8 ±6.1 


8.040 ±0.440 


29.0 ±2.9 


2.05 ±0.82 


H (50.0) 


37.5 ± 


0.6 


-6.82 ±0.95 


0.0894 ±0.0059 


47.8 ±8.9 


14.478 ±0.859 


42.3 ±5.3 


0.93 ±0.28 



Series I: q = 0.1, Af+ = 20,000 



model (fi) 


<*cc> 






(rc) 


(a/;) 


(«min> 


(A/"esc) 




A (1.25) 


3859.0 ±00.0 


-8.62 ±0.00 


0.0021 ±0.0000 


33.0 ±0.0 


1.000 ±0.000 


422.0 ±0.0 


0.52 ±0.00 


B (1.5) 


3318.0 ± 


0.0 


-9.09 ±0.00 


0.0030 ±0.0000 


36.0 ±0.0 


1.016 ±0.000 


292.0 ±0.0 


0.79 ±0.00 


C (2.0) 


2236.0 ± 


0.0 


-5.64 ±0.00 


0.0044 ±0.0000 


23.0 ±0.0 


1.081 ±0.000 


140.0 ±0.0 


0.34 ±0.00 


D (3.0) 


1094.0 ± 


0.0 


-3.84 ±0.00 


0.0114 ±0.0000 


34.0 ±0.0 


1.170 ±0.000 


55.0 ±0.0 


0.56 ±0.00 


E (5.0) 


608.0 ± 


0.0 


-4.21 ±0.00 


0.0103 ±0.0000 


16.0 ±0.0 


1.770 ±0.000 


65.0 ±0.0 


2.82 ±0.00 


F (10.0) 


213.0 ± 


0.0 


-3.70 ±0.00 


0.0241 ±0.0000 


24.0 ±0.0 


2.950 ±0.000 


38.0 ±0.0 


1.34 ±0.00 


G (25.0) 


78.0 ± 


0.0 


-4.95 ±0.00 


0.0863 ±0.0000 


88.0 ±0.0 


7.680 ±0.000 


30.0 ±0.0 


1.41 ±0.00 


H (50.0) 


48.5 ± 


0.0 


-7.92 ±0.00 


0.0447 ±0.0000 


25.0 ±0.0 


14.046 ±0.000 


48.0 ±0.0 


0.78 ±0.00 



Series II: q = 0.1, A/"* = 2500, INS 



model (fj.) 


(tec) 




<*min> 


(rc) 




(A/"esc) 


(E CS c ) 


A (1.25) 


623.9 ±11.1 


-5.64 ±0.10 


0.0098 ±0.0005 


15.2 ±0.5 


57.5 ±2.8 


0.89 ±0.22 


B (1.5) 


501.5 ± 


9.1 


-5.05 ±0.10 


0.0133 ±0.0008 


14.4 ±0.8 


42.0 ±2.5 


1.30 ±0.56 


C (2.0) 


304.4 ± 


8.0 


-4.09 ±0.12 


0.0225 ±0.0018 


15.2 ±1.0 


16.9 ±1.1 


0.35 ±0.18 


D (3.0) 


139.4 ± 


4.6 


-3.51 ±0.07 


0.0330 ±0.0023 


13.5 ±0.9 


8.5 ±0.9 


0.15 ±0.11 


E (5.0) 


73.6 ± 


3.0 


-3.51 ±0.08 


0.0437 ±0.0026 


12.3 ±0.8 


8.4 ±0.8 


1.43 ±0.40 


F (10.0) 


32.4 ± 


1.7 


-3.58 ±0.13 


0.0722 ±0.0059 


16.3 ±1.9 


5.9 ±0.8 


1.01 ±0.29 


G (25.0) 


18.3 ± 


0.9 


-5.28 ±0.33 


0.0845 ±0.0070 


16.3 ±2.9 


13.6 ±1.5 


0.78 ±0.11 


H (50.0) 


15.2 ± 


0.9 


-5.52 ±0.27 


0.0892 ±0.0050 


13.9 ±1.4 


25.7 ±2.2 


0.82 ±0.10 



Series II: q = 0.1, A/"* = 2500, OUT 



model (^t) 



(tec) 



<*min> 



(rc) 



(€min> (A4sc) 



(E CS c) 



A (1.25) 
B (1.5) 
C (2.0) 
D (3.0) 
E (5.0) 
F (10.0) 
G (25.0) 
H (50.0) 



673.7 ± 8.6 
646.1 ±11.3 
508.5 ±10.8 



304.1 ± 
153.7 ± 
92.7 ± 

52.5 ± 

47.6 ± 



7.3 
1.7 
3.0 
1.5 
2.3 



-6.08 ±0.16 
-5.12 ±0.14 
-4.16 ±0.12 
-3.23 ±0.05 
-2.98 ±0.06 
-3.23 ±0.08 
-4.58 ±0.27 
-6.86 ±0.29 



0.0088 
0.0129 
0.0227 
0.0476 
0.0756 
0.0892 
0.0857 
0.0664 



±0.0004 
±0.0008 
±0.0014 
±0.0026 
±0.0042 
±0.0033 
±0.0037 
±0.0032 



14.8 ±0.6 

17.1 ±0.9 

18.2 ±1.2 
24.8 ±1.8 
29.6 ±2.8 
23.8 ±1.6 
17.5 ±1.2 

9.7 ±0.7 



61.4 ±2.5 
50.1 ±2.5 
32.6 ±1.6 
13.0 ±0.8 
8.5 ±0.6 
10.0 ±0.9 
17.0 ±1.6 
32.0 ±2.8 



1.03 ±0.41 
0.66 ±0.14 
1.77 ±0.42 
1.10 ±0.30 
1.40 ±0.32 
2.17 ±0.47 
1.76 ±0.26 
1.94 ±0.32 



Series III: 


q = 0.05, A/"* 


= 2500 












model (fi) 


(tec) 


(*min) 


(rc) 




<£min> 


(A^csc) 


(Eesc) 


K (1.25) 


689.3 ± 7.8 


-5.88 ±0.21 


0.0094 ±0.0005 


16.0 ±0.6 


1.000 ±0.000 


60.7 ±2.5 


0.53 ±0.12 


L (1.5) 


605.3 ±11.8 


-5.19 ±0.19 


0.0140 ±0.0009 


17.2 ±0.8 


1.001 ±0.000 


51.0 ±2.7 


0.61 ±0.16 


M (2.0) 


428.8 ± 9.8 


-3.89 ±0.11 


0.0264 ±0.0022 


19.2 ±1.4 


1.009 ±0.008 


25.8 ±1.7 


0.82 ±0.28 


N (3.0) 


219.3 ± 7.4 


-3.20 ±0.06 


0.0443 ±0.0028 


19.5 ±1.6 


1.092 ±0.034 


9.3 ±1.1 


0.85 ±0.27 


P (5.0) 


96.1 ± 4.5 


-3.07 ±0.06 


0.0672 ±0.0046 


23.6 ±2.4 


1.510 ±0.034 


4.3 ±0.7 


1.60 ±0.54 


Q (10.0) 


44.9 ± 1.4 


-3.34 ±0.11 


0.0896 ±0.0053 


23.5 ±2.4 


2.638 ±0.068 


4.2 ±0.6 


1.06 ±0.19 


R (25.0) 


27.6 ± 1.2 


-4.53 ±0.21 


0.0903 ±0.0046 


19.0 ±1.9 


4.348 ±0.290 


9.0 ±1.4 


1.37 ±0.44 
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Series IV: q = 0.2, jV* = 2500 



model (fj.) 


(ice) 




<<t>mi n > 


(rc> 




<5min) 


<A/" C sc) 


(Eesc ) 


T (1.25) 


657.1 ±11.0 


-5.88 ±0.17 


0.0102 ±0.0006 


15.2 ±0.9 


1.001 ±0.000 


62.5 ± 2.6 


0.81 ±0.21 


U (1.5) 


549.4 ± 


7.6 


-5.00 ±0.14 


0.0134 ±0.0009 


14.6 ±0.8 


1.000 ±0.002 


49.5 ± 2.7 


1.17 ±0.28 


V (2.0) 


343.8 ± 


9.5 


-4.51 ±0.11 


0.0166 ±0.0011 


12.3 ±0.6 


1.067 ±0.012 


26.7 ± 2.2 


1.04 ±0.35 


W (3.0) 


174.0 ± 


3.9 


-3.87 ±0.11 


0.0244 ±0.0020 


11.8 ±0.8 


1.334 ±0.016 


16.9 ± 1.3 


0.61 ±0.04 


X (5.0) 


86.3 ± 


2.2 


-3.55 ±0.09 


0.0391 ±0.0028 


12.0 ±0.9 


2.069 ±0.059 


13.9 ± 0.8 


0.81 ±0.08 


Y (10.0) 


48.6 ± 


2.4 


-3.97 ±0.16 


0.0564 ±0.0041 


11.9 ±1.0 


4.063 ±0.161 


21.1 ± 1.6 


1.06 ±0.11 


Z (25.0) 


22.7 ± 


1.0 


-4.72 ±0.24 


0.0869 ±0.0049 


15.5 ±1.6 


9.676 ±0.267 


30.1 ±1.97 


0.74 ±0.07 



Series V: q = 0.4, A/"* = 2500 



model (fi) 


(tec) 




<4-mi„> 




(AQ 


«min> 




(E csc ) 


Tf (1.25) 


639.7 ± 


9.9 


-5.79 ±0.16 


0.0092 ±0.0007 


14.0 ±0.8 


1.000 ±0.000 


63.0 ±2.4 


0.68 ±0.13 


U/ (1.5) 


523.2 ± 


8.2 


-5.24 ±0.12 


0.0117 ±0.0006 


14.4 ±0.8 


1.007 ±0.003 


51.4 ±2.7 


0.81 ±0.22 


V/ (2.0) 


372.3 ±10.3 


-4.95 ±0.15 


0.0149 ±0.0007 


13.3 ±0.6 


1.094 ±0.013 


42.3 ±2.9 


0.85 ±0.24 


W/ (3.0) 


212.8 ± 


1.3 


-4.27 ±0.10 


0.0195 ±0.0014 


11.3 ±0.6 


1.436 ±0.039 


38.7 ±2.4 


0.57 ±0.03 


X/ (5.0) 


130.7 ± 


1.0 


-4.19 ±0.14 


0.0274 ±0.0015 


10.4 ±0.6 


2.200 ±0.047 


43.7 ±2.9 


0.81 ±0.05 


Y/ (10.0) 


66.8 ± 


2.5 


-4.26 ±0.12 


0.0440 ±0.0026 


9.9 ±0.5 


4.316 ±0.087 


50.9 ±3.5 


0.88 ±0.05 
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